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presented in the final paper. 


INTRODUCTION 

The motivation of this effort is based on two primary elements j. First, 
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Ica^ simui^ion^l'rbe of^ significant value toward design. This will help to improve the 
safety of future space missions. 
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speed up will be due^o enhancements in computer hardware platforms. Remaining por 



tion of the speed-up will be contributed by advances in algorithms and by efficient parallel 
implementations. The following section outlines the initial effort and steps taken in order 
to meet this milestone. 



Figure 1. Schematic view of RSTS boost pump and its components. 


REUSABLE LAUNCH VEHICLE (RLV) TURBOPUMP INDUCER 

Rotational Speed : 7850 RPM 
Mass Flow: 9093 G PM 
Re : 7.99e+7 



Geometry Surface Pressure 


Figure 2. Geometry and surface pressure for the boost pump inducer 


APPROACH AND RESULTS 

The current geometry for the LOX boost turbopump has various rotating and stationary 
components, such as inducer, stators, kicker, hydrolic turbine, where the flow is exteremly 
unsteady. The components of the entire turbopump are plotted in Figure 1 . Figure 2 
shows the geometry and computed surface pressure of the inducer. The inducer and the 


hydraulic turbine rotate in different rotational speed. This causes severe unsteady interac- 
tions between rotating and stationary parts. To handle the geometric complexity and mov- 
ing boundary problems, overset grid scheme is incorporated with the solver that new 
connectivity data is obtained at each time step. The Chimera overlapped grid scheme 
allows sub-domains move relative to each other, and provides a great flexibility when the 
boundary movement creates large displacements. In order to validate the current 
unsteady solution procedure, SSME shuttle upgrade pump configuration has been 
selected for validation purpose. Figure 3 shows the geometry of the test rig for this pump 
being tested at NASA-MSFC facilities. In this particular configuration, SSME impeller is 
unshrouded. 



Figure 3. Geometry of SSME shuttle upgrade pump impeller 

The development of the MPI, hybrid MPI/Open MP versions of the INS3D code is com- 
pleted, and MLP version of the code is currently underway. The serial version of INS3D 
code is a multidimensional incompressible Navier-Stokes solver based on overset grid 
technology. The algorithm is based on the pseudocompressibility method as developed 
by Chorin 1 . The pseudocompressibility algorithm introduces a time-derivative of the pres- 
sure term into the continuity equation; the elliptic-parabolic type partial differential equa- 
tions are transformed into the hyperbolic-parabolic type. The original version of the INS3D 

code 2 with pseudocompressibility approach utilized the Beam-Warming approximate fac- 
torization algorithm 3 and central differencing of the convective terms. Since the convec- 
tive terms of the resulting equations are hyperbolic, upwind differencing can be applied to 
these terms. The current versions of the INS3D code 4 use flux-difference splitting based 
on the method of Roe 5 . The upwind differencing leads to a more diagonal dominant sys- 
tem than does central differencing and does not require a user specified artificial dissipa- 
tion. The viscous flux derivatives are computed by using central differencing. In the 
steady-state formulation, the time derivatives are differenced using the Euler backward 
formula. The equations are solved iteratively in pseudo-time until the solution converges 
to a steady state. In the time-accurate formulation, the time derivatives in the momentum 


equations are differenced using a second-order, three-point, backward-difference for- 
mula. After the discretization in time, the pseudocompressibilty term and pseudo-time 
level are introduced to equations. The equations are iterated to convergence in pseudo- 
time for each physical time step until the divergence of the velocity field has been reduced 
below a specified tolerance value. The matrix equation is solved iteratively by using a 

nonfactored Gauss-Seidel type line-relaxation scheme 6 , or GMRES scheme 7 which 
maintains stability and allows a large pseudo-time step to be taken. Details of the numer- 
ical method can be found in Reference 4. 
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Figure 4. Computed surface pressure for SSME-HPFT impeller. 


INS3D-MPI 8 is based on the explicit massage-passing interface across processors and is 
primarily suited for distributed memory systems. The primary strategy is to distribute the 
zones across a set of processors. During the iteration, all the processors would exchange 
boundary condition data between processors whose zones shared interfaces with zones 
on other processors. A simple master-worker architecture was selected because it is rela- 
tively simple to implement and it is a common architecture for parallel CFD applications. 
All I/O was performed by master MPI process and data was distributed to the workers. 
After the initialization phase is complete, the program begins its main iteration loop. Com- 
puted surface pressure of the shrouded SSME impeller is shown in Figure 4. The solution 
is obtained by using INS3D-MPI with line-relaxation scheme. Computational grid for this 
case has 24 zones and 2.8 Million grid points. The MPI version of the code has coarse 
grain parallelism, and the number of processors is limited by the number of zones. Figure 
5 shows Mflops for this computation on SGI Origin 2000 platform. The average speedup, 
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OpenMP hybrid version of the code OpenMP processors from hybrid parallel 
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Figure 5 INS3D-MPI performance in the SGI Q2K, SSME-HPFT Impeller. 



Figure 6 1NS3D-MPI /Open-MP performance in the SGI 02K, SSME-HPFT Impeller. 

The third parallel version of the c ??® ^ technique develop^ NASA Ames 
shared memory Multi-Level Paral ' e ’ ^i into INS3D code. INS3D-MPL code has inter 
Research Center has been in ^ or P°J a ^ ' n '° ' „ e lism is obtained by using OpenMP d.rec- 
zone parallelism with MLP, and '^ra zon P a computations has been completed 

fives Computational model for the overset grid surfaces, 

by using an overset grid technology. Figure f snows 




This particular impeller con,i 9“^ ,ion ,u a s . ™ outehonaTm^ zones wi,h 35 

inlet guide vane arij^ffuse^Wades, the P inte , ide vanes and diffuser blades. 

Million grid points. Impeller blades roraxe code Qn 512 proce ssors 

This configuration will be co^puted wit Ce nter. Since it is an unsteady simu- 

SGI Origin 2000 platform at NASA Am ®\7®^ C !J oraae fea ture extraction, and overset 
lation, some of the additional algori £ m developed at NASA-Ames is 

connectivity synchromzatio . P . (OVERFLOW-D 11 ) will be run 

incorporated in the procedure. Overset connective code (OVERFLOW , 

simultaneously with the flow solver. 



Figure 7. Computation grid of SSME impeller with 60 zones 


and 19.5 Million grid points. 
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Figure 8. Flowchart from CAD to flow feature extraction tor the current procedure. 

Figure 8 shows a typical flowchart tor the for the entire tur- 
"S S -as wher/the advancements are currency under- 

way are outlined. 
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